Low-mj^ 7r+-7r Asymmetry Enhancement from Hadronization of QGP 
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We show that in sudden hadronization of QGP a non-equilibrium value of the pion phase space 
occupancy parameter 7, > 1 is expected in order to accommodate the entropy excess in QGP and the 
process of gluon fragmentation. When 7, is near to its maximum allowed value, pion overabundance 
is shown to arise at low mx, where the charged pion asymmetry is also amplified. These effects 
are considered and we show that their magnitude suffices to explain pertinent experimental data 
obtained in 160-200A GeV Pb- and S- induced reactions. 



Highly relativistic collisions of heavy nuclei are un- 
der experimental and theoretical study of which the 
primary aim is the discovery of the deconfined quark- 
gluon plasma (QGP) phase and the understanding 
of its transformation into hadrons. We shall demon- 
strate how the measurement of low central rapid- 
ity abundance of charged pion yields, dN{'K~)/dm± and 
dN{'K~^)/d'mj_, provides a sensitive way to probe the 
outcome of hadronization of QGP. We will, in particu- 
lar, show that the relative deformation of low-m^ pion 
spectra provides a measure for chemical (particle) non- 
equilibrium overabundance phenomenon. 

Thermal equilibration (i.e., equilibration of the energy 
distribution) of final state hadrons is seen in particle spec- 
tra to be a well working hypothesis. This, however, does 
not imply chemical (i.e., abundance) equilibration. Since 
chemical hadron abundance equilibrium is approached on 
the same time scale as the duration of the heavy ion col- 
lisions, it has been proposed to analyze the abundances 
of final state hadrons allowing chemical non-equilibrium 
phase space occupancies ||^. This has been applied ini- 
tially only to describe the abundance of strange quarks 

, which were early on understood to be not necessarily 
produced rapidly enough ||]. More recently, it has been 
suggested that chemical equilibration of light quarks and 
gluons is indeed not assured at RHIC energies and be- 
yond |^-|]. On the other hand, at SPS, the excess of 
pion multiplicity has been observed and understood to 
result from excess entropy content of the primordial high 
entropy phase formed in the interaction 0,^. This leads 
naturally to the hypothesis of light quark overabundance 
generated, e.g., by gluon fragmentation in hadronization 
of the QGP high entropy phase. This can be most easily 
described by counting light valence quark abundance in 
hadrons ]s|-[ll|], as was done for strangeness. 

Indeed, when the multi-particle production processes, 
in 158-200A GeV S- and Pb-Pb coUisions carried out 
at CERN-SPS, have been analyzed using the methods of 



the statistical Fermi model |T^, light quark overabun- 
dance and the associated chemical non-equilibrium was 
found. Overall, the results of this analysis show several 
features that combined do make a strong case that at 
least the strange hadronic particles seen at CERN-SPS, 
are emerging nearly directly from sudden hadronization 
of deconfined QGP phase of hadronic matter |11|. More- 
over, the recognition that the hadron abundances are 
better described with oversaturation of the quark phase 
space allows a notable reduction of the chemical freeze- 
out temperature, Tf = 145 ± 5 MeV, accompanied by 
a reduction of the baryochemical potential, fis = 200- 
210 MeV, is consistent with the underlying dynamical 
picture. What remains open is the issue if the predom- 
inant hadronic fraction, pions, also follow the sudden 
hadronization model, a point we address here. 

This general discussion introduced already a few sta- 
tistical parameters of the Fermi model of well established 
meaning and we will complete the list now. Aside of the 
chemical freeze-out temperature Tf, further chemical pa- 
rameters control hadron abundances: we note the light 
quark fugacity and and the phase space occupancy 
of light quarks 7^ — 7„ = 7^. We recall that 7^ controls 
overall abundance of quark pairs, while Ai control the 
difference between quarks and anti-quarks of given fla- 
vor. The difference between A^ and 7^ is that, for quarks 
and anti-quarks the same factor 7^ applies, while the an- 
tiparticle fugacity is inverse of the particle fugacity. The 
phase space occupancy 7^ is determined to exceed unity 
significantly, once introduced: results obtained switch- 
ing on progressively more parameters of the Fermi model 
show considerable improvement in the capability of the 
model to describe the data |p^ . 

The proper statistical physics foundation of 7^ is ob- 
tained considering the maximum entropy principle: it has 
been determined that while the limit 7^ — > 1 maximizes 
the specific chemical entropy, this maximum is not very 
pronounced, indicating that a system with dynamically 



1 



evolving volume will, in general, find more eff'ective paths 
to increase entropy, than offered by the establishment of 
the absolute chemical equilibrium ]T^ . 

It is usual to introduce the geometric average of the 
light quark fugacities, related to the baryochemical po- 
tential: 



(1) 



Similarly, one introduces the chemical potentials for up 
and down quarks: 

/irf = T In Arf , ^„ = Tin A„ , S^i = fid ~ Mu • (2) 

The difference in the number of net up (u — u) and net 
down (d — d) quarks participating in the dense matter 
fireball is given by the initial condition, i.e., which nu- 
clei collide. In a hadron gas, 5fi/ fj,q varies as function of 
chemical freeze-out temperature as shown in Fig. 1] . 

The relative number of particles is controlled by fu- 
gacity and phase space occupancies of the constituents. 
The composite hadronic particle chemical occupancy and 
fugacity is expressed by constituent contributions: 
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(3) 



The abundances of particles produced in QGP disinte- 
gration is most conveniently described by considering the 
phase space distribution of particles: the abundance of 
a hadron h, with j-components, freezing out from the 
source is 



(4) 



For example, the effective chemical fugacities for tt^ 
composed of light quark- antiquark pair of a differing fla- 
vor is 



77r- = 7, 



I Ad 
1 1T' ' 



ln+ = Iq 



2 



Ad' 



while for neutral pions we simply have: 



1^° = In ■ 
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The case of pions is, however, exceptional in another 
way: the SPS data analysis shows that the pion yield is 
governed by a large fugacity, and thus it is imperative 
to revert in the Fermi model to use Bose distribution 
function: 

(Pn^ 1 
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= ^ml+p^. (7) 



We see that the range of values for 7^ is bounded from 
above by the Bose singularity: 



In < % 



When 7^ the lowest energy state (in the contin- 

uum limit with p ^ ) will acquire macroscopic occupa- 
tion and a pion condensate is formed. Such a condensate 
'consumes' energy without consuming entropy of the pri- 
mordial high entropy QGP phase j^,^ and we cannot 
expect this to occur in hadronization of a high entropy 
phase into a low entropy confined phase. However, a 
sudden hadronization process will naturally have the ten- 
dency to approach the limiting value 77r 7^ in order to 
more efficiently connect the deconfined and the confined 
phases, since, as we show in Fig. ^ the entropy density 
is nearly twice as high at — 7^ than at 7,^ = 1- We 
note that all thermodynamic quantities of interest here 
do not develop a singularity in the limit 7,^ ^ 7^. 
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FIG. 1. Dependence of pion gas properties A'^/V-particle, 
E/V-eneTgy and S/V-entTopy density, as function of 7, (for 
Ad = A„) at T = 142 MeV. 

We thus consider as a possible and indeed likely QGP 
hadronization mechanism, the conversion of the excess 
of entropy of the deconfined state into oversaturation of 
pion phase space. This mechanism can replace the more 
conventional approaches that also described the 'absorp- 
tion' of the excess QGP entropy: 

1) formation of a mixed phase which allows the volume of 
the fireball of hadronic matter to grow in the hadroniza- 
tion; the excess of QGP in entropy is converted into the 
greater volume of lower entropy density hadronic gas. 

2) reheating which allows the entropy excess to be accom- 
modated in the momentum 'volume' of the phase space. 

Of course, both these effects can coexist, and also they 
can, in principle, coexist with our proposed third mech- 
anism, the chemical non-equilibrium. In fact, which ex- 
act path is taken is only found in a microscopic descrip- 
tion which allows for all these effects with appropriate 
relative rates. A model study of the dynamics of the 
hadronization process shows that the reaction mechanism 
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here proposed is possible |15|, though a full understand- 
ing of the microscopic hadron phase transformation re- 
mains at present beyond our capability. However, aside 
of the fact that the data analysis heavily favors our pro- 
posal, we note that the conventional processes are by far 
more complex and thus more difhcult to realize on the 
short time scale 0(10"^^ s) available in heavy ion col- 
lisions, as compared to a rather straight forward gluon 
fragmentation, accompanied by recombination of quarks 
into pions, which is the probable microscopic origin of 
the large chemical nonequilibrium value of 7^. 

Let us return now to the study of pion abundances: the 
slight chemical asymmetry originating in the larger neu- 
tron number than proton number in heavy nuclei, leads 
to light quark abundance asymmetry, Nd > iV„, and re- 
quires that one of the charged pion fugacities (tt") in- 
creases (towards the condensation point), while the other 
charged pion (tt"*") fugacity is moved further away from it. 
The result of this is clearly visible in Eq. (|^) when com- 
bined with Eqs. (||,^): as function of energy this slight 
change has the effect of enhancing the yield of low energy 
7r~ over that of tt^ . There is a much more significant ef- 
fect when we are near to 7^ than it is the case for 7^ = 1 
which was explored earlier p6| . In the following, we will 
address in detail the statistical parameters as they are 
determined for the 158~200AGeV interactions and pre- 
sented in table ^. We note that for the Pb-Pb system, 
the value of = 1.10 compensates the Coulomb effect 
distortion of the phase space of strange quarks in QGP. 
Thus, it is a measure of the magnitude of the Coulomb 
potential at hadronization, and hence of the freeze-out 
volume [|lO| . We so estimate the volume of the fireball at 
chemical freeze-out to be F ~ 2100 fm'^, and this value is 
used to check if the absolute yields of pions, and in par- 
ticular here negative hadrons (7r~,K~,p), are consistent 
with measurement of the experiment NA49 ||l^ . We find 
within our approach that among the 200 'negatives' at 
central rapidity Ay = 1, there are (roughly) about 145 
directly produced vr^ , 20 p and K~ , and only about 30 
iT~ that arise from resonance decays. 

Fig. H presents a case study of the effect on the thermal 
pion TO_L spectra. We integrate the longitudinal momenta 
in the range of rapidity —0.5 < y < 0.5 and consider the 
difference in charged pion yield divided by the sum. 
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as function of m±, for small m± — m^r < 0.3 GeV, show- 
ing the result on logarithmic scale. The lowest line shows 
the small effect of isospin asymmetry in Pb-Pb collisions 
when the pion phase space is governed by the equilibrium 
chemical condition 7^ = 1; this result corresponds, for the 
range of our statistical parameters here applicable, to the 
results reported in Rcf. jl^. As we move up, for each of 
the lines in Fig.|^, we increment 7g = 1 by 0.1 and the 
thick and top line corresponds to the value implied by 
the analysis of the experimental results, 7^ — 1.6. The 



magnitude of the effect is now quite surprising (remem- 
ber logarithmic scale), and most importantly, we do see 
a very clear sensitivity to the value of 7^. While the 
'modulator' of this sensitivity is the isospin asymmetry 
Xd/^m this quantity is in principle measured by the re- 
sults shown in Fig.^, since for m± — m„ > 0.3 GeV we 
have: 
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It must be stressed that the result shown in Fig. Q must 
be compared to experiment with some caution, since it is 
just a case study aiming to show the extraordinary sen- 
sitivity of the D'^(rn±) observable to the quantity 7^. In 
order to be able to interpret actual data, one has to study 
the diluting pion component from resonance decay and 
the effect of collective flow. It can be hoped that the flow 
effect is small, and indeed cancelling in the pion yield ra- 
tio. The contribution of resonance decays turns out also 
to be rather small in the range of statistical parameters 
considered in view of the small chemical freeze-out tem- 
perature. 
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FIG. 2. Thermal pion asymmetry m±_ spectrum, for 
7g = 1 to 1.6 (thick line) by step of 0.1, integrated in the 
rapidity window —0.5 < y < 0.5. 



TABLE I. Chemical freeze-out statistical parameters for 
Pb-Pb 158^ GeV data ^ and S-Au/W/Pb 200^ GeV data 
[^. Last line is adapted from j^. Fig. 1]. 

Pb S 



Tf [MeV] 
\ 

\s 

7s/7<7 



142 ± 3 
1.61 ± 0.02 
1.10 ± 0.02 

1.6 ± 0.5 
0.80 ± 0.05 
0.07 



144 ± 2 
1.51 ± 0.02 
1.00 ± 0.02 
1.41 ± 0.08 
0.69 ± 0.03 
0.04 



3 



Given the magnitude of the 7g-effect, we also studied 
in Fig.^ the spectrum of neutral pions, which is practi- 
cally indistinguishable from the average of the charged 
pion spectra, (7r+ +7r~)/2. The solid line shown in Fig.^ 
are not renormalized, and (on logarithmic scale) the dif- 
ference between the curves for m± — m^r > 0.3 GeV is 
the result of the varying yield coefhcient 7^. More im- 
portantly, for m_L — m,r < 0.2 GeV with growing 7^ we 
see a strong up-bending of the pion spectrum. Dashed 
lines in Fig. || show results for relatively small Ay = 0.02 
and demonstrate that there is just a minimal spectral 
deformation due to a finite central Ay = 1 integral. 
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FIG. 3. 



Thermal tt spectrum. 
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7q = 1 to 1.6 by step of 0.1, integrated in the rapidity window: 
-0.5 < y < 0.5 (full lines), and -0.01 < y < 0.01 (dashed 
lines), the latter normalized to the result with Ay = 1 at 
m± — rrin = 0.5 GeV. 

In order to understand the experimental pion spectra 
it is necessary to consider the 7g-effect, combined with 
the flow effect and the hadron disintegration feeding of 
pion spectra. We recall here that aside of direct influ- 
ence on primordial pions, 7^ also impacts the relative 
yield of primordial mesons compared to the yield of sec- 
ondary mesons from resonance decays. Thus omission 
to consider 7q-effect can mislead the data interpretation. 
We recall here a recent analysis of the 7r° spectra which 
has reached somewhat unorthodox conclusions about the 
freeze-out conditions as seen in Ref. Tablet]. We be- 
lieve that allowance for chemical non-equilibrium, specif- 
ically the introduction of 7^ ~ 7^ should alter the con- 
clusions of the WA98 collaboration. Even though they 
study the pion spectrum at rather large m_L, 7^ impacts 
the cocktail mix of direct with decay pions. There is an- 
other detailed fit which also restricted the single particle 
pion spectra to be exactly at chemical equilibrium [p^ . 
However, these authors do note that in order to describe 
the pion yield, they need a pion fugacity 7^*^" = 1.36, 



a value which may (barely) allow their tacit assumption 
7g(T — 100 MeV) = 1, made when they interpret par- 
ticle spectra. At the higher value of chemical freeze-out 
temperature a range of 7^ is expected that in our opin- 
ion requires reevaluation of their results. On the other 
hand, the parallel analysis of Schlei and collaborators 
pO| , which did not attempt a precise fit, and did not 
describe the low m± pion spectra, and which obtains a 
similar chemical freeze-out temperature as used here, ap- 
pears to us to remain valid. 
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FIG. 4. Central rapidity charged pion ratio as function 
of m± : data are by NA44 collaboration j2l] , arbitrarily nor- 
malized to unity at high m±. The lines are our result for 
7g = 1.6, also renormalized to unity at large m±. 

It would be interesting to see what effect would arise 
after explicit allowance for a pion fugacity is made in 
these studies. Since we do not model here the effect of 
collective flow as was done in Refs. [^-|o|, we cannot 
ourselves answer this question in full. However, there is 
data from NA44 on charged pion ratios ||2l[]. Considering: 

1) the relatively small dilution of primordial pion spectra 
by pions from resonance decays, expected for the range 
of chemical freeze-out parameters shown in table ^ and 

2) the cancellation of the flow deformation of the spectra 
in spectral ratios of particles of the same mass, 

we have evaluated the expected jq-eSect for the charged 
pion ratio results presented by NA44 collaboration pi[ , 
as shown in Fig.^. We note that we slightly overpredict 
at small m± these results, an effect that is required given 
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the inclusion of hadron resonance decay feeding into the 
experimental data. 

Several attempts have been made before to explain 
these NA44 data. The experimental group presented 
calculational results obtained with RQMD(vl.08) 
which could not account for the data. Even though NA44 
considered in their numerical model also hadron reso- 
nance decay effect, a scheme has been proposed to create 
the observed pion asymmetry from strangeness asymme- 
try ||2^ : in presence of finite baryon number there are 
more s-carrier baryons than antibaryons and more s- 
carrier mesons than s-carrier mesons. Of course, such 
a description can easily be falsified by study of strange 
hadron abundances. In thermal models, for 7^ = 1, 
the NA44 asymmetry cannot be explained for the widest 
range of reasonable statistical parameters, a fact we have 
realized long ago, and which is implicit in other work 
p6t . These studies have thus lead to the believe that 
Coulomb field effects in single particle spectra are the 
source of pion asymmetry and a large number of authors 
has explored this hypothesis |2^-|2^]. It is quite possi- 
ble that in part the effect of the pion asymmetry is due 
to the Coulomb effect, and in part explained by the 7^- 
effect here described, considering that the value 7^ = 1.6 
obtained from the study of hadron abundances has a sig- 
nificant statistical error. 

The Coulomb effect work is based on the picture of 
charged pions acted upon by the Coulomb field after for- 
mation at a classically well determined point. The pion 
asymmetry is a result of the hypothesis that the pion for- 
mation weight is actually equal for both 7r+ and tt^ 
at equal local momentum: 
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(11) 



where i?„ is a static or dynamic formation surface, 
the WKB approximation: 



In 
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We believe that considerable improvement in these stud- 
ies of the Coulomb effect is possible. In our view, 
the Coulomb field influences the quantum density of 
states obtainable from the scattering phase r/. This in- 
troduces in addition to the statistical density matrix 
weight Tr e^^^ an additional density of states weights 
drf ~/dE, drf^/dE, which due to the presence of attrac- 
tive (for 7r~) and repulsive (for tt^) Coulomb potential is 
Coulomb deformed at asymptotic pion momenta values 
that sample the physical size of the potential. We note 
that for Rn — 8-10 fm this effect is expected to occur for 
— TO^ < tt'^ /2Rf^mTt ~ 10-20 MeV in static case. In 
quantum formulation the presence of (Coulomb) poten- 
tial alters locally the reduced wave length X = ?i/p^(a;), 
but not the asymptotic energy E^^ of a quantum state. 
The work most closely related to this approach [ p3|j26| 
takes Coulomb wave-functions (thus better than the 



WKB illustration made above) but adds ad hoc density 
of state hypothesis, see Ref Eq. (2)], using here the 
density of matter in configuration space folded with the 
Coulomb wave. The pion asymmetry is arising as result 
of substitution of Eq. ([l2|) in Eq. (^. 

Our findings have interesting bearing on the forthcom- 
ing experimental RHIC results. The data from the ex- 
periment PHOBOS |2^ should in our opinion also allow 
to observe the 7q-effect, both in soft pion spectra, and in 
pion asymmetry, provided that the non-equilibrium oc- 
cupancy, 7g ~ 1.6, which we associated with the QGP 
hadronization, also occurs at RHIC — as we expect [ p8| . 
The analysis of these data should follow the pattern out- 
lined here: from the slight inequality of the abundance at 
moderate m± , according to Eq. , the value of chem- 
ical asymmetry can be deduced and this is used to de- 
termine the value of 7^, given the shape of the low- rn^ 
spectral asymmetry following from Eq. (^. We note that 
Coulomb _analysis also predicts a similar effect at RHIC 
energies 
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We have shown that the phase space overoccupancy 
7g ~ 1.6 of pions arising in the sudden hadronization 
process of the deconfined quark-gluon plasma enhances 
the low-rnj_ pion spectra, and amplifies the abundance 
ratio of charged pions. We have quantitatively demon- 
strated the great sensitivity of these effects to the value 
of 7g > 1 and have shown how this 7q -effect helps under- 
stand SPS soft pion results of NA44 and WA98 collab- 
orations. We conclude that the reported excess of soft 
pions and the asymmetry of soft charged pion spectra 
provides evidence for sudden and explosive hadroniza- 
tion (glue fragmentation) of the entropy rich deconfined 
phase. This finding about the mechanism of hadroniza- 
tion corroborates the conclusions we obtained studying 
high m± strange hadron spectra. 
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